function y = genhumps(x)
% GENHUMPS  function
% x0 = [506.2,506.2,..., 506.2]
n=length(x);
y=0;
for i=1:n-1;
    y=y+sin((2*x(i))^2)*sin((2*x(i+1))^2)+0.05*((x(i))^2+(x(i+1))^2);
end


x0=1:10;
y0=1:10;
[x,y]=meshgrid(x0,y0);
z=sin((2*x)^2)*sin((2*y)^2)+0.05*(x.^2+y.^2);
mesh(x,y,z)


